In-Medium Effects on Charmonium Production in Heavy-Ion Collisions 
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Charmonium production in heavy-ion collisions is investigated within a kinetic theory framework 
incorporating in-medium properties of open- and hidden-charm states in line with recent QCD lattice 
calculations. A continuously decreasing open-charm threshold across the phase boundary of hadronic 
and quark-gluon matter is found to have important implications for the equilibrium abundance of 
charmonium states. The survival of J/tp resonance states above the transition temperature enables 
their recreation also in the Quark-Gluon Plasma. Including effects of chemical and thermal off- 
equilibrium, we compare our model results to available experimental data at CERN-SPS and BNL- 
RHIC energies. In particular, earlier found discrepancies in the ip' /ip ratio can be resolved. 



The production systematics of heavy-quark bound 
states in (ultra-) relativistic collisions of heavy nuclei (A- 
A) is believed to encode valuable information on the hot 
and dense strong-interaction matter formed in these re- 
actions . Based on the notion that charm quark pairs 
(cc) are exclusively created in primordial (hard) nucleon- 
nucleon (N-N) collisions, it has been suggested [2j that a 
suppression of observed J / ip mesons in sufficiently central 
and/or energetic A- A reactions signals the formation of 
a deconfined medium (Quark-Gluon Plasma=QGP), as 
tightly bound cc states are conceivably robust in hadronic 
matter. While theoretical Q and (indirect) experimen- 
tal evidence jj] supports N-N collision-scaling of total 
charm production, it has recently been realized H, IE 
that coalescence of c and c quarks can induce significant 
regeneration of charmonium states in later stages of A- 
A collisions, especially if several pairs are present (e.g., 
iV ce =10-20 in central Au-Au at RHIC). This is a direct 
consequence of (elastic) c-quark reinteractions, facilitat- 
ing the backward direction of charmonium dissociation 
reactions, J/ip + X x ^ X 2 + c + c (D + D). 

Recent lattice computations of Quantum Chromody- 
namics (QCD) have revealed important information on 
charm(onium) properties at finite temperature T, most 
notably: (i) an in- medium reduction of the open-charm 
threshold which is surprisingly continuous even across the 
phase transition region ||, and, (ii) the survival of JJjb 
and r\ c states as resonances in the QGP phase [jj lioj . 
with essentially unmodified masses. Note that the use of 
charmonium spectral functions, in connection with ap- 
propriate asymptotic states, incorporates both (static) 
screening and (dynamical) dissociation mechanisms. 

In this letter, we propose an approach that implements 
charm properties inferred from lattice QCD as micro- 
scopic in-medium effects into a kinetic rate equation. It 
enables a comprehensive treatment of charmonium disso- 
ciation and regeneration across the phase transition, con- 
necting hadronic and QGP phases in a continuous way, 



not present in previous calculations. E.g., in the kinetic 
approach of Ref. 0, J ftp's were only considered in the 
QGP phase using their vacuum binding energy, whereas 
recent transport calculations 0, do not invoke the 
notion of a (equilibrium) phase transition. In addition, 
in-medium open-charm masses will significantly affect 
the equilibrium levels of charmonium, as employed in sta- 
tistical models 0, • We thus conceptually improve our 
earlier constructed "two-component" model |l3| . which 
combined statistical production at hadronization with 
suppression in hadronic and QGP phases. Employing 
a schematic thermal fireball model, the rate equations 
are applied to heavy-ion reactions at SPS and RHIC. 

We begin by evaluating in-medium masses of charm 
states. Charmonium masses appear to be essentially 
unaffected at finite T , even above T c (critical tempera- 
ture) HE3, which we assume from now on. Unquenched 
QCD lattice studies of the free energy, Fqq(t,T), of a 
heavy-quark pair however, exhibit a plateau value 
for r — > oo which gradually decreases with increasing 
T. In the hadron gas (HG), this can be interpreted as 
a reduction of D-meson masses (supported by QCD sum 
rules 0|) being driven by a reduced constituent light- 
quark mass, m*, due to (partial) chiral symmetry restora- 
tion ^^|. We implement this also for other open-charm 
hadrons by employing a Nambu-Jona-Lasinio model cal- 
culation for m* at finite T and quark chemical poten- 
tial \i q ^(| 5 neglecting heavy- light quark interaction and 
kinetic energies. For small fi q , the typical light-quark 
mass reduction amounts to Am 9 (T c )~140 MeV. To as- 
sess open-charm states in the QGP, a key element is the 
essentially continuous behavior of the open-charm thresh- 
old through T c . Within a quark description above T c , we 
propose to identify this with an in-medium charm-quark 
mass, m*~1.6-1.7 GeV [HI- In fact, m*(T c ) is fixed such 
that the total open-charm density smoothly matches the 
hadronic side. The difference to the bare charm-quark 
mass (m c ~1.3 GeV) is naturally attributed to a thermal 
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correlation energy of heavy quarks in the QGP. 

Next, we determine the thermal widths (or inverse life- 
times), r^ = (r^,) _1 , of charmonium states tp. They are 
related to (inelastic) dissociation cross sections, cr^ ss , by 
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where i runs over all matter constituents with ther- 
mal distributions Within the QGP, parton- 
induced "quasifree" break-up, i + tp — > i + c + c 
(i=g,q,q), is utilized 0. For small binding energies, 
ER=2m*-m,i, <Aoc n , the latter (with a s ~0.25) has been 
shown 01 to dominate standard gluo-dissociation (lif . 
g + tp — > c+ c. In addition, essentially unbound charmo- 
nium resonances such as t/>' and x-states can be treated 
on an equal footing. The thermal J/tp width turns out to 
be ~100 MeV at T~250 MeV, reminiscent to (quenched) 
lattice results ^3 under comparable conditions. 

For J ftp's in the HG, we compute inelastic cross sec- 
tions with 7r- and p-mesons within a flavor- > S'/7(4) 
effective lagrangian formalism [20ll2l| . including in- medi- 
um charm hadron masses [2^ as described above. The 
free x an d i 3 ' cross sections are approximated by geo- 
metric scaling [23J, supported by quark-model calcula- 
tions [2^[. If the DD threshold moves below the char- 
monium mass, we additionally evaluate direct decays, 
tp — > Z?l), accounting for wavefunction-node effects ac- 
cording to Ref. . Despite the increase in available 
phase space due to reduced masses, the charmonium life- 
times in the HG remain longer than in the QGP. 

In-medium charm properties have important conse- 
quences for charmonia in heavy-ion reactions. To illus- 
trate this point, consider the J/tp number in thermal 



equilibrium, 
density 



■Vfb 1 ^ (Vfb- 3-volume), with the 
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Here, we have included chemical off-equilibrium effects 
via the charm-quark fugacity j c which is adjusted to the 
total number, N cBl of cc pairs in the system via 0,0 



1 h{l c N op , 
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V FB ^2 n^(T,j c ) (3) 
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where N op =VFBn op (m* c D \T) denotes the total 
equilibrium-number of open-charm states in the thermal 
fireball. This procedure resides on the expectation that 
essentially all charm-quarks are created in primordial N- 
N collisions, i.e., N c5 does not chemically equilibrate 0. 
Thermal equilibration of charm quarks, however, is 
conceivable |27| . implying a hierarchy of relaxation 



times as r^g 



> T FB > T% e 



The resulting J/tp 

equilibrium abundances, N^ q , are shown in Fig. ^ under 
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FIG. 1: Equilibrium J/tp abundances (with m^,(T) = ni^ ac ) 
in an isotropic, adiabatically expanding system at fixed N C c 
for SPS and RHIC conditions for two values of the charm- 
quark mass above T c , and for free and in- medium charmed 
hadron masses below T c . 



conditions resembling central Pb-Pb (Au-Au) collisions 
at SPS (RHIC). One observes a large sensitivity to the 
(in-medium) open-charm masses. At fixed N cEl larger 
values for to* (or m* D ) imply a thermal suppression 
of open-charm states so that an increasing number 
of anti-/charm quarks is redistributed into charmonia 
(cc states). Enforcing continuity on the open-charm 
spectrum across T c then has the interesting consequence 
that, due to the volume increase in the hadronic phase 
(reducing j c ), the equilibrium charmonium level on the 
QGP side is significantly larger than on the HG side [l7j . 
i.e., J/tp formation is favored in the QGP. 

Also note that, due to a constant N c5 , the equilibrium 
J/tp numbers in the hadronic phase grow with decreas- 
ing temperature. Thus, if J/tp's equilibrate close to T c , 
subsequent hadronic reactions will tend to increase their 
abundance (which has indeed been found in transport 
calculations at RHIC energies ^10]), quite contrary to 
the commonly assumed hadronic dissociation. 

To model the time dependence, AT^(r), of charmonia 
in heavy-ion collisions, we utilize the dissociation (forma- 
tion) time t^, Eq. and equilibrium density, Eq. J3J), 
within a kinetic rate equation, 



dNjP 
dr 
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The underlying temperature and volume evolution, T(r) 
and Vfb(t), are taken from a (schematic) thermal fire- 
ball expansion which is consistent with observed hadro- 
chemistry and radial flow characteristics (as employed 
before for thermal dilepton production H^)- 

The simple form of the gain term in Eq. (0J resides 
on the assumption that the surrounding light and open- 
charm constituents (quarks or hadrons) are in thermal 
equilibrium. Following Ref. |23| . we relax this assump- 
tion for c-quarks by introducing a thermal relaxation 
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time correction, 1Z = 1 — exp( 



duces Nf? in the early phases 



■/ dr/r^ 6 ™), which re- 
^ i±i uik ctU ij Varying T * herm within a 

factor of 2, the thermal yield is affected by ±10%. 

Another correction concerns the effective volume over 
which the charm quantum number is conserved (figuring 
into Eq. via the Bessel functions). Clearly if only few 
cc pairs are present, their pointlike primordial production 
implies that they cannot explore the entire fireball vol- 
ume in the early stages j^. This problem is well-known 
from strange particle production at fixed target energies, 
where a phenomenological "correlation volume", Vo, has 
been introduced to localize strangeness conservation [i^ . 
We adopt the same procedure here for local charm con- 
servation by replacing Vfb{t) in the argument of the 
Bessel functions in Eq. © with Vb(Y)=47r(ro + (v c )t) 3 /3. 
ro—1.2 fm represents a minimal radius characterizing the 
range of strong interactions, and (i; c )~0.5c is the aver- 
age relative speed of produced c and c quarks as inferred 
from experimental D-meson p^-distributions [30] . We 
checked that our results are not very sensitive to the pa- 
rameterization for Vo (e.g., varying (v c ) within a factor 
of 2 requires a change in the effective coupling constant 
a s by ±30%, to fit SPS data). 

The rate Eq. (0J is readily integrated over the fireball 
evolution of the collision for J/ip, ip' and x' s i once ini- 
tial conditions are specified. In an A- A reaction at given 
centrality, N c5 is determined by the collision-scaled N- 
N cross section. The initial charmonium numbers are 
also assumed to follow hard production (i.e., empirical 
fractions of N cB observed in N-N collisions), subject to 
(pre-equilibrium) nuclear absorption with a recently up- 
dated absorption cross section <7 a 6 s =4.4 mb jUl- Typical 
values for the fireball thermalization time range from i 
(RHIC-200) to 1 fm/c (SPS). 

We first confront our approach to NA50/NA38 data 
in ^5^=17.3 GeV Pb-Pb at CERN-SPS. Fig. H dis- 
plays the ratio of J/tp^fifi to Drell-Yan dimuons as a 
function of centrality. The agreement between model 
(solid line) and data is fair for semi-/central collisions 
(at Dr>100 GeV, the data can be reproduced by ac- 
counting for transverse energy fluctuations and losses in 
the minimum bias analysis at impact parameters close 
to zero [3^, cf. Ref. [23 )■ Since the initial J/ip num- 
ber is well above the equilibrium level (cf. lower arrow in 
Fig. nj, J/ip regeneration (the gain term in Eq. J3J) is 
very moderate (dot-dashed line). Therefore, in line with 
our previous findings J/ip suppression is the main 
effect at SPS energies. 

In peripheral collisions, the suppression appears to be 
slightly overestimated. We believe this discrepancy to 
reside in limitations of our fireball description. In par- 
ticular, thermalization is expected to be delayed (and/or 
incomplete) at large impact parameters due to less en- 
ergetic initial conditions. This is also borne out of hy- 
drodynamic models, which, e.g., reproduce the observed 
elliptic flow for mid-central collisions, but overestimate it 
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FIG. 2: Centrality dependence of J/?/>/Drell-Yan dimuons 
at SPS; NA50 data |3l| are compared to our results with 
(solid line) and without (dashed line) charmonium regenera- 
tion (represented by the dot-dashed line) . The dotted line in- 
cludes longer thermalization times to in peripheral collisions. 
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FIG. 3: Centrality dependence of the ip' /ip data |34{ at SPS 
compared to our results with (full line) and without (dashed 
line) in-medium reduced D-meson masses. 



for peripheral ones. A suitable increase of the equilibra- 
tion time by up to a factor of ~3 |33| indeed improves 
the agreement at small Et, cf. dotted line in Fig. El 

Another important observable is the ip'/tp ratio. In 
Ref. [2^], the ip' dissociation rates were too small by a 
factor of ^5 to account for NA50 data ^34]. With in- 
medium D-meson masses, however, T^f d increases sub- 
stantially, primarily due to the opening of the ip' — > DD 
decay channel. As a result, the ip' /ip data are reasonably 
well described, cf. Fig. |21 

We finally examine the impact of in-medium modi- 
fications at RHIC, cf. Fig. Since reduced D-meson 
masses entail a lower J ftp equilibrium level (cf. Fig. 0, 
the regeneration of J/tp's is somewhat less pronounced 
than the statistical production with free hadron masses 
in the 2-component model [U [3f|. Nevertheless, in 
central Au-Au collisions, regenerated J/ip's still exceed 
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FIG. 4: J/ip yield per binary N-N collision versus partici- 
pant number in y / ijv r iv=200 GeV Au-Au collisions. Prelim- 
inary PHENIX data ,36j are compared to our model calcu- 
lations; dot-dashed line: thermal regeneration; dashed line: 
suppressed primordial production; band: total J ftp yields 
with different values for in-medium open-charm masses. 



the suppressed primordial contribution, with the total 
yield (solid curve) in line with first PHENIX data [3r| . 
The uncertainty in the charm-hadron mass reduction 
is illustrated by the band in Fig. corresponding to 
80 MeV<Am g (T c )<250 MeV with accordingly matched 
charm-quark masses in the QGP. In any case, most of 
the regeneration occurs above T c . 

To summarize, we have proposed a conceptually im- 
proved model of charmonium production in heavy-ion 
collisions. In-medium modifications of both open-charm 
and charmonium states have been modeled in accordance 
with recent finite-T QCD lattice calculations. The ap- 
parent reduction of the open-charm threshold with in- 
creasing T has been linked to (partial) chiral symme- 
try restoration in the SU (2)-sector via decreasing con- 
stituent light-quark masses in charmed hadrons. The 
continuity of the open-charm threshold across the phase 
transition has been encoded in thermal charm-quark 
masses in the QGP phase, together with T-independent 
charmonium masses and J/ip resonance states surviving 
above T c . Upon inclusion of chemical and thermal off- 
equilibrium effects, we solved kinetic rate equations for 
the time evolution of charmonia in heavy-ion reactions. 
The results for J/ tp centrality dependences compare well 
with data from SPS and RHIC. The earlier predicted 
transition from suppressed production at SPS to pre- 
dominant (thermal) regeneration at RHIC persists, with 
both mechanisms residing on QGP formation. A much 
improved description of the iji'/ip ratio at SPS emerged 
due to hadronic in-medium effects, accelerating tp' disso- 
ciation. Forthcoming measurements by NA60 (SPS) and 
PHENIX (RHIC) will be of great importance to scruti- 
nize the proposed approach. 
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